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, We study the generation of non-Gaussianity in models of hybrid inflation with two inflaton fields, 

C^l ■ (2-brid inflation). We analyse the region in the parameter and the initial condition space where a 

(■^ ' large non-Gaussianity may be generated during slow-roll inflation which is generally characterised 

, by a large /nl, tnl and a small gNL- For certain parameter values we can satisfy tnl 3> / NL . The 

■ bispectrum is of the local type but may have a significant scale dependence. We show that the loop 

corrections to the power spectrum and bispectrum are suppressed during inflation, if one assume 
that the fields follow a classical background trajectory. We also include the effect of the waterfall 
field, which can lead to a significant change in the observables after the waterfall field is destabilised, 
depending on the couplings between the waterfall and inflaton fields. 

O !■ PACS numbers: 98.80.Cq 
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I. INTRODUCTION 

The increasingly accurate observations of the cosmic microwave background [l[ motivate the study of inflation 
beyond linear order in perturbation theory. Since there are many models of inflation which give similar predictions 
for the spectral index of the scalar perturbations and the tensor-to-scalar ratio 0, Q it is important to consider higher 
, order observables. Non-Gaussianity has emerged as a powerful discriminant between different models of inflation. 
■ Any detection of primordial non-Gaussianity would rule out the simplest models of single field inflation. 

Many ways to generate a large non-Gaussianity have been proposed in the literature 0, H, @, 0, IS HI • If the inflaton 
has a non-canonical kinetic term then the inflaton field perturbations may already be non-Gaussian at Hubble exit 
qq ' during inflation [lol ]. Otherwise the perturbations are Gaussian at Hubble exit and any non-Gaussianity must then be 
generated on super Hubble scales [ill, Q ■ Popular methods to generate a large non-Gaussianity after inflation include 
having a feature in the inflaton potentia l fl3ll. t he curvaton scenario [Til ], modulated reheating/preheating TEI. 16j, an 
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inhomogeneous end of inflation [T7|, lis llSl, l20j, |2l| or by loop domination of the bispectrum or trispectrumj22L ]2c 
Work has also been done on calculating the non-Gaussianity during inflation with a separable potential [2~il. |25|. l26l. l27l| 
and with more general potentials [28l . [23 | . Recently the authors of this paper have shown that it is possible to generate 
- . i a large non-Gaussianity during multiple field inflation while keeping all of the slow-roll parameters much less than 
unity [30| . 

In this paper we study this possibility of generating a large non-Gaussianity during slow-roll inflation in detail for the 
specific model of hybrid inflation. We consider carefully the parameter constraints and initial conditions required to 
do this. Unlike in the previous work [3(| we here also consider the effect of the waterfall field required to end inflation. 
We can therefore also compare and contrast our work to a recent paper on generating a large non-Gaussianity at the 
end of inflation [2l[. Depending on the values of the couplings between the two inflaton fields and the waterfall field 
observable quantities may change when the waterfall field is destabilised. This change in the primordial curvature 
perturbation is possible because there are isocurvature perturbations present at the end of inflation. 
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We also calculate further observables for this model, such as the trispcctrum and the scale dependence of the 
bispectrum. In most cases the trispectrum is large through a large tjvl whenever the bispectrum is large but the 
relation between them also depends on the initial conditions. 

Recently it has also been shown that for two-field hybrid inflation it is possible to generate a large loop correction to 
the bispectrum and trispectrum [22l.l23j. For a special trajectory it was shown that the bispectrum can be observably 
large even when the tree level value of the bispectrum is slow-roll suppressed. Here we show that this is not possible 
for a classical trajectory, since we require that the classical motion of the field dominates over its quantum fluctuations 
in order that a slow-roll calculation is valid. 

The plan of our paper is as follows: In Section [TT] we review the generation of a large bispectrum during hybrid 
inflation and calculate the trispectrum. In Section IIIII we include the waterfall field and consider its effect at the end 
of inflation for different couplings between the waterfall and inflaton fields and further evolution after inflation. In 
Section IIVI we calculate the scale dependence of the non-Gaussianity parameter. We conclude in Section [V] In the 
Appendix, we show that the loop correction is subdominant to the tree level term for the power spectrum, bispectrum 
and trispectrum. 



II. NON-GAUSSIANITY DURING HYBRID INFLATION 

We consider a model of two field hybrid inflation, whose potential is given by 



which is vacuum dominated, i.e. which satisfies \r\ w ^ 2 \ <C Mp and |% X X 2 | *C Mp. We assume that inflation ends 
abruptly by a waterfall field which is heavy during inflation and hence doesn't affect the dynamics during inflation. 
In this section we calculate observables during slow-roll inflation. We will consider the full potential including the 
waterfall field in Sec. IIIII We will see that this can lead to a change in observables on the surface where the waterfall 
field is destabilised. 

In the vacuum dominated regime the slow-roll solutions are 

<p = tp*e-*>* N , X = x , e -^ N , (2) 

where '*' denotes the value at the horizon exit. Throughout this section whenever we write a quantity without making 
it explicit at which time it should be evaluated, we mean the equation to be valid at any time N e-foldings after 
Hubble exit and while slow roll is valid. Generally we will be interested in quantities at the end of inflation, in which 
case we take N = 60. 

The slow-roll parameters are 

_ 1 2 V> 2 _1 2 _ , ,a\ 

We note that the dominant slow-roll parameters r] vv and r\ xx are constants during inflation in the vacuum dominated 
regime and that they are much larg er than the slow-roll parameters e v and e x throughout inflation. 

Using the 5N formalism [HI, l32l. l33l. l34l . f35j ] (for details see the appendix) we can calculate the power spectrum 
Vq, spectral index tiq, tensor-to scalar ratio r, where Vt is the power spectrum of tensor perturbations, and the 
non-linearity parameter /nl in this model pH [24. l25l. [30L |36| . at leading order, 
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According to [3(j, large non-Gaussianity can be realised in either of two regions 

cos 2 ^^^-— ^— <1, or sin 2 0=— ^— ~^^<1. (8) 

</> + X £ v + e x V + X <v + e x 

Since both regions are symmetrical [30| (before specifying the values of 7/ w and ?y xx ), in the rest of this paper we will 
focus on the second region (Region B or D in Ref. [3(|). In this region where e v e x , |/nl| > 1 is fulfilled by the 
condition, 



sin 2 ^ <sin 4 4\/^^-l h (9) 
6sin 2 6> J K ' 

in other words, 

IVxxl' 1 ^ 4 ^^ £ sm 2 9~ ^ < \Vxx\- (10) 
This condition implies three conditions on the parameter 9: 

• 2„* 1 f5\ 2 {3V, ,2 . 2n 5, sin 2 6> 24 1 

sin 2 0*<- - - \Vxx\ , sin 2 < - \ri YX , — 5 — > — r . 11 

3V6/ W 6 1 /xx 1 ' sin 2 6>* 5 |?7 XX | V ' 

Note that in this region sin^ ~ VxxXI 1 {tIvp^P) ■> f rom Eq. ([!]) we require N(r] vv — rj xx ) > 1 so that sin 2 9 grows 
significantly during inflation. 



A. Simplified formula for the observables when /nl is large 

We can substantially simplify all of the above formula in the case where /nl is large. We define the quantity 

f= (™Y I (™Y = € JL e ^ vv -r, xx )N (12) 

\dx*J \dLp*J e v 

In the region we are considering where /nl is large, this is approximately given by the initial and final angles of the 
background trajectory with different exponents 



sin 4 (9 
sin 2 9* 



(13) 



We note that f can be either larger or smaller than one. 
In the case of large non-Gaussianity it follows that 



v ^2^A 1+ t/ {n)N ) = r\^ ■ (14) 



(15) 
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The first condition in Eq. pip implies that 



— <1. (18) 



We therefore require a very small value of x* m order to have a large non-Gaussianity. We will see in sec III Bl that 
there is a limit to how small we can make x* while x still follows a classical slow-roll trajectory. In practise for 
relatively large values of rj vv — rj xx this places a bound on how large we can make /nl- 

The sign of /nl is determined by the sign of r] xx . The amplitude of /nl depends exponentially on the difference of 
the slow-roll parameters, V(f>(p~ r lxx' which we require to be positive to be in the branch of large non-Gaussianity where 
sin 2 6 <C 1. However the spectral index depends on the weighted sum of the slow-roll parameters, so it is possible to 
have a large non-Gaussianity and a scale invariant spectrum. However it is not possible to have a large and positive 
/nl and a red spectrum of perturbations. We will see in Sec. IIIIBI that by including the effect of the waterfall held 
this conclusion may change, depending on the values of the coupling constants between the two inflaton fields and 
the waterfall field. 

Using the SN formalism we can also calculate the non-linearity parameters which parameterise the trispectrum. 
There are two shape independent parameters, which may be observationally distinguishable and are given by [37l ] 



5 N AB N A N B N AB N BC N A Nc 25 N ABC N A N B N C 

6 (NcNcf ' ™ L ~ {NnNnf ' 9NL ~ 54 { N D N D f ; 



where summation over the fields is implied over the repeated indices A, B, C, D. For comparison we also give the 
formula for /nl [35j [- In the regime where \gNb\ > 1 it is given by 

IQ ^fapy -27fr x ) -?7 XX 

9NL = — 7—1 /NL- (2U) 

3 1 + r 

Details of the calculation for both g^L and tml are in the appendix, in particular see Eqs. (fSTj) . (|52"|) and (|67p . We 
see that g^L is subdominant to /nl and hence won't provide a competitive observational signature. The current 
observational bound on the local type of the bispectrum from WMAP data is — 9 < /nl < HI at the 2a level [l[ 
and currently there is no observational constraint on the trispectrum. If there is no detection it is expected that with 
Planck data the bounds will be reduced to about |/nl| ^ 10 and tnl < 560 [3^ ]. 

However in the regime where tn l is large, which is similar to the region where /nl is large we find 

From Eq. (|17|) we see that for given model parameters fy w and 7; xx the value of /nl is maximised for f = 1, i.e. when 
we choose the initial field values to satisfy 



V* Vxx 

In this case the non-linearity parameters are given by 



(22) 



/NL = ^ xxe 2W( ^- Ito) , r^ = 2(j/ NL ) . (23) 

However tjvl is maximised at a slightly different point, for f = 1/2. 

It follows from (f2"Tj) that ttvl > (6/nl/5) 2 , so ttvl m ay be large and provide an extra observable parameter for 
this model. This inequality between tjvl and /nl is true in general [l6j . and equality is reached whenever a single 
field direction during inflation generates the primordial curvature perturbation. However it is usually assumed that 
t nl ~ /nl since both arise from second derivatives in the SN formalism. In fact for our model it is possible to have 
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TABLE I: Table showing some initial conditions for the hybrid inflation model that lead to large levels of non-Gaussianity. 
The end point of fields, f , the bispectrum and trispectrum non-linearity parameters, spectral index and tensor-to-scalar ratio 
are shown. They are evaluated when the number of e-foldings from the end of inflation is Nk = 60. 



a small /nl (and hence also a small Qnl) but a large and potentially observable t^l- For this we require that f <C 1, 
in practice if we make it too small it may no longer be possible to satisfy the classical constraint (|2~4"|) discussed in 
Sec. Ill Bl In the final example in Table [J we give an explicit example of parameter values which give rise to an /nl 
which is probably too small to be detected with Planck but with a very large trispectrum through t/vx > 10 3 that 
should be detectable at a high significance. For another example with Jnl,9nl iS 0(1) but tnl 1, see [39[. In 
contrast it has been shown in several papers [53, S3, EH, 0, [43| that in the curvaton scenario where the curvaton 
has a non-quadratic potential it is possible to realise \gNh\ 3> 1 while tnl — (6/nl/5) 2 is small with some tuning 
of parameters. Recently it has also been shown that in a single field model with a non-canonical kinetic term it is 
possible for certain tuned parameters to achieve a much larger trispectrum than bispectrum [44j . We note that the 
non-Gaussianity in their model has a very different shape dependence to the k independent non-linearity parameters 
we are considering in this paper. 

In Table HI we give some explicit examples of values of ?7 W , tj xx , ip* and x* which lead to a large non-Gaussianity. 
Using Eq. ([5]) we also calculate the spectral index. The contours of /nl of potential for a specific choice of -q vv and 
r\ xx is given in Fig. 5 in our previous paper [3(|. We correct tensor-to-scalar ratio in Table 1 of [30| . 1 The first 
example in the Table [J shows that it is possible to have |/nl| — 100 and a scale invariant spectrum. We also see that 
it is possible to generate a large non-Gaussianity during slow roll with r] vv and rj xx both positive or both negative, 
or when one is positive and the other negative corresponding to a saddle point. 



B. Can the loop correction dominate? 

Cogollo et al. have calculated the effect of the loop correction to the primordial power spectrum and bispectrum (23 
and even more recently Rodriguez and Valenzuela- Toledo have calculated the loop correction to the trispectrum [23j . 
In both cases this was for the special case of a straight background trajectory where one of the fields is zero throughout 
inflation. This loop correction arises from taking into account the contribution to the power spectrum and bispectrum 
arising from terms in the SN expansion which are non- leading in the expansion of the field perturbation Sip*, see 
e.g. [43. However they can still be significant if the coefficient to the term given by the second derivative of N with 
respect to the subdominant field is extremely large and the leading term for the same field is small or zero, e.g. [46| . 
These "higher order" terms are usually neglected, but has shown the first explicit example of an inflation model 
where they cannot be neglected. They consider a 2-field hybrid model with the same potential as Eq. ([T]) in the special 
case of an unstable straight trajectory along one of the axes, with r\ Vip and rj xx both negative. In this case, they find 
(for certain initial values) that one of the loop correction terms is dominant over the tree level term and can generate 
an observable /nl- 

However if the value becomes too small then the motion of the \ field will become dominated by quantum fluc- 
tuations rather than the classical drift down the potential, 3Hx — ~~W X = —Wor) xx x/Mp, which we have assumed. 
In order that we can neglect the effect of the quantum fluctuations, the condition we require on the background 
trajectory is [13] 

\x\ir/H 2 > V3?2. (24) 



1 The value of the tensor-to-scalar ratio in Table 1 of |30| is incorrect. The correct value for the first row is 0.005, for the second row is 
0.026 (as given in Table [I} and for the third row r = 0.002. This does not change the discussion or conclusions of |30| . 



6 



We require the condition above to be satisfied throughout inflation. Since the fields in our model either increase or 
decrease monotonically it is sufficient to check that the constraint is satisfied both at Hubble exit and at the end of 
inflation. In the appendix we show that satisfying this condition requires that the loop corrections are suppressed. 
We have checked that for the examples in Table Q] this condition is satisfied. We stress that this suppression of the 
loop corrections follows from the requirement that the background trajectory of the inflation fields is dominated by 
the classical motion. However it would be interesting to investigate the large loop corrections found in [H, using a 
calculation which includes the quantum corrections to the x field, which is set to zero throughout inflation in [22l . [23| . 

Because we have shown that the loop correction is always suppressed compared to the tree level terms we are 
justified in using the formula for the tree- level non-Gaussianity parameters as given by Eq. (|19|) . 



III. NON-GAUSSIANITY AFTER HYBRID INFLATION 



In this section we include the effects of the waterfall field p which is required to end hybrid inflation. Inflation ends 
when the waterfall field is destabilised, i.e. when its effective mass becomes negative. During inflation the waterfall 
field is heavy and it is trapped with a vacuum expectation value of zero, so we can neglect it during inflation. The 
end of inflation occurs when the effective mass of the waterfall field is zero, which occurs on a hypersurface defined 
in general by [13, HI, 

a 2 = G(<p, X ) = SiV + g\x\ (25) 
which is realised by the potential W((p, x), defined by Eq. (p}, where Wo is given by 

W^\G{^X)P 2 + \{p 2 ' ■ (26) 

Here g\ (32) is the coupling between the ip (x) field and the waterfall field. In general the hypersurface defined by 
this end condition is not a surface of uniform energy density. Because the SN formalism requires one to integrate up 
to a surface of uniform energy density we need to add a small correction term to the amount of expansion up to the 
surface where the waterfall field is destabilised, which we will consider in sec. IIII CI 

We note here that the hybrid potential, which we have written in the form of a sum potential can also be written 
as a product potential in the limit of vacuum domination. We show this so that we can use the formula from [2lj | . we 
consider a potential of the form 

W(<p tX ) = W exp (^Vw-Jtfs) ex P (^*Xj§jr) ' ( 2? ) 

We note that the slow-roll parameters are not exactly the same from these two ways of writing the potential, but that 
we are only interested in calculating results at leading order in slow roll, and in this limit the two ways of writing the 
hybrid potential are equivalent. 

This is an example of a model with an inhomogeneous end of inflation, i.e. where inflation ends at slightly different 
times in different places. This has been studied as a method for converting the inflationary isocurvature perturbation 
into the primordial curvature perturbation [l§| . It has also been shown for the hybrid potential we are considering 
that this can be used to generate a large amount of non-Gaussianity, for certain parameters values and fine tuning of 
the parameters [III [2l|. However these papers concern the large non-Gaussianity generated at the end of inflation 
rather than during slow- roll inflation, by having a very large ratio of couplings 31/(72 <C 1. Here we consider the case 
where g\ and 32 have the same order of magnitude with g\j g\ = Vvv/Vxx m sec. IIII Al and with g\ = g\ in sec. IIII Bl 



A. gl/gl = rivf/Vxx 

In this case we have chosen the coupling constants (which can satisfy g\ < and/or g\ < 0) such that the surface 
where the waterfall field is destabilised corresponds to a surface of uniform energy density. In this case the value of all 
observable quantities such as the power spectrum and non-Gaussianity are the same as those we calculated previously 
which were valid at the final hypersurface of uniform energy density during inflation. 

With these values of the coupling constants, large non-Gaussianity can only be generated with the special conditions 
we have outlined earlier in this paper and it is given by Eq. (|14[) - Eq. Q17p . We have checked that our result for 
/nl, Eq. Q in this case is consistent with the formula for /nl in [2l| . This provides a check on the algebra. We 
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FIG. 1: Left: The contour plot of /nl which have the end condition same as uniform energy density hypersurface g\jg\ = 
Vifif/Vxx in IHI Al in the plane of phi and chi, which denote the values of the fields when the given scale leaves the horizon. The 
values of /nl, 0, 10 and 20, are shown on the corresponding contour. For example with the given end point, the trajectory is 
shown as a line denoted by "Traj". The cross point with Nk = 60 line is the value of phi and chi when the scale corresponding 
to Nk = 60 leaves horizon. In this scale the value of /nl is around 9. Right: Directly from the left figure, we can read the 
scale dependence of /nl on the given trajectory. For the given trajectory the spectrum of /nl is shown in the figure. For this 
trajectory we used the third example in Table |H i.e. r) vv — 0.08, rj xx = 0.01, and the end point is fixed so that the field values 
<p* = 1 and x* = 0.0018 lead to the number of e-foldings Nk = 60. 



emphasise that this large non-Gaussianity is not from the end of inflation but from the evolution during inflation. To 
compare between our paper and [2l[ we note that in their notation m\ = ?y w , m| = r\ xx , = — m\m2<p e XeW / (2Z) 
and e x = mim\ip e x e Y / (2X) where the parameters W,X,Y and Z are defined in [2l| . 

In the left hand plot of figureQ] we show the contour plot of /nl in the (if, x) plane of the field values when the given 
scale leaves the horizon, /nl increase as the number of e-folding from the end of inflation increases, which is shown 
in the right figure. This can be easily understood from Eq. (| 1 T[) . For a small c-folding number, Nk, (corresponding to 
small scales) f is also small and /nl is small. As Nk increases, f increases as well as /nl- However for f > 1, which 
corresponds to Nk > [\og(e v /e x )]/2(ii vv — T] xx ) for a given end point, /nl approaches its the maximum value, which 
corresponds to 

f^ X = l € fv xx - (28) 
u t x 

In the example given in Fig.[T] this happens around Nk > 120, thus the maximum is not shown in the right-hand plot 
of Fig. [1] in the given range of Nk- The value of Nk which gives the maximum value can be modified by the condition 
of end of inflation and can be approached Nk ~ 60 which is shown in the next section and Fig. [2] 



B. gl = gl 

In this case, the end of inflation given by the condition in Eq. (|25[) does not occur on a uniform energy density 
hypersurface [2l| . In this subsection we will show how the non-Gaussianity is modified by the condition at the end of 
inflation with this specific example. In general we expect there to be some modification to non-Gaussianity from the 
end of inflation, except in the special case we considered in IIII Al 

In this case the power spectrum and the non-linearity parameters are (2l| 

v Wo (^K v )e~ 2 ^ N + je x h\ x )e-^ N 

c 247r 2 Mj, (ev/Vw + ex/Vxx) 2 
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(32) 



16e* 



1 + 



(33) 



r nl 



(34) 



We can see that the observables at the end of inflation are changed by the additional ratio of iJ^/Vxx before f. 
Furthermore, /nl has an additional factor of rj w /rj xx . 

In Table [ill we show the values of /nl, uq — 1 and r for the same parameter values as we used in Table H The first 
two examples show that if r) vv = — r\ xx then the observables are unchanged except that the sign of /nl is switched. 
The second example in the table shows that in this case it is possible to have a red spectral index and a positive value 
of /nl- For many values of the initial parameters when \r\ w jr\ xx \ ^ 1 the magnitude of /nl decreases compared to 
Tablcfl] We see that this is the case for most examples given in the table. This is most strikingly illustrated in the fifth 
example in the two tables. We can see by comparing (JT7J) and (|3"4")l that in the case where f < 1 and \il<p<p/ r }xx\ ^ 1 
that the magnitude of /nl will decrease by about \Vipv/' l lxx\ 3 ■ However for most values of f and the ratio TJfp/Vxx 
the change is much smaller. In some cases /nl increases, by the largest amount when f < 1 and \r]<p<p / Vxx\ > 1- 
example of this where |/nl| grows by more than an order of magnitude is shown in the sixth column of Tables [J and 
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TABLE II: Same as Table [T] but with different end condition, g\ = g\ as used in sec. IIII Bl 



Due to the additional factor before r, /nl is increased at small scales compared to the case in sec. IIII Al and decreased 
at large scales due to the additional factor. Thus the maximum /nl is smaller than that of sec. IIII Al by Ti vv /rj xx . At 
small scales (where f < 1), /nla = (v vv /Vxx)f^B- At l ar g e scales (where f > 1), / NLA = (% x /?7 w )/nlb- 



C. The end of hybrid inflation 

As we mentioned in the previous sections, the end of inflation given by the condition in Eq. |25]) does not occur on 
a uniform energy density hypersurface (2lT |. i.e. the energy density is slightly different for the different end points. For 
the SN formalism to be valid, SN should be calculated at the uniform energy density hypersurface. This is explained 
in (2^1 with the extra term N c in Eq. (3.14) of [2(j. They assume that the universe has become radiation-dominated 
right after inflation. This means that at the time t e when the water fall field starts to become unstable, the energy 
density of the inflaton fields changes to radiation instantly conserving the energy density, ignoring the details of the 
water fall field dynamics. If we define the time t c sometime in the radiation dominated era, the change in the number 
of e-foldings between t e and t c is 



N c = - In 
4 



Wj_ 

W 

This can be understood like this: During radiation domination, [4 

dp 



(35) 



which gives 



m -4,, (36) 



N c = N(t c )-N(t e )= 1 -ln^= 1 -ln^l= 1 -ln^l+ 1 -ln^. (37) 
v ' y ' 4 Pc 4 W c 4 W 4 W c y ' 



where we can put W c = Wo which anyway does not change SN. This is Eq. (3.15) in [201 . We require that the 
perturbation from this term, SN C , is smaller than from 5Nf. This was shown for the hybrid potential with a linear 
exponential dependence on the field in [2(| , but it was not shown in [2l[ where they consider the same hybrid potential 
as we do, with a quadratic dependence on the fields. Therefore we need to check the validity of neglecting <5iV c in our 
two examples in the previous subsections. 

In the first example with g\j g\ = rj vv / l r\ xx , as in sec. IIII Al the hypersurface at the end of inflation is a uniform 
energy density hypersurface, which means SN C = 0. Thus in this case, there is no further change of Jnl, if we assume 
an instantaneous change to radiation. 

In the other case with g\ = g\ , as in sec. IIII Bl wc obtain 



Nc = \ ln ^k = \ {jlvV ^ !+ ^xX 2 ) ■ 08) 



Taking the perturbation of N c , we find 



Wc - TrziVxx -rhpv) [(sin2 7 )(<Si7 + ^7) + (cos 2 7 )(<5i 7 ) 2 ] , (39) 
°9i 



where 7 is defined by tan7 = x/f a ^ the end of inflation satisfying Eq. (|25|) . Using the expressions for <5i 7 and ^27 
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((Eq (A2) and (A4) respectively in [2l|), this becomes 

— 77^— + W,n<n— sin 9^i/ Vxx ( ~77T~ 



sin 27 



?7 XX tan7 + 77^/ tan 7 2 77 xx tan7 + tan 7 

x2l ( 4 °) 

1 { (-'Vxx/ cos2 7 + »W shl2 7) sin 2 7 + cos 2 7 (?7 XX tan 7 + 77 w / tan 7) } [ri xx 5 -^f - r\ w ^ 

2 (77 xx tan7 + 77 V¥ ,/tan7) 3 

The perturbation of N from the flat hypersurface at horizon exit during inflation to the hypersurface at the end of 
the inflation (with a corrected factor of one half missing in (A. 8) in [2lj), 

l/tan 7 ^+tan 7 ^ t -1/ tan 7 (^) ' + tan 7 (^) ' l/(sm 7 cos7) (y xx ^ - (41) 
77 xx tan7 + 77 w /tan 7 2 7y xx tan 7 + 77 w / tan 7 (t? xx tan 7 + 77 w / tan 7 ) 3 

The final SNf from the flat hypersurface at horizon exit during inflation to the uniform energy density hypersurface 
during radiation dominated era is the sum of both 

SN f = 5N + SN C . (42) 

By comparing the corresponding term of Sip*/ip, S\*/x e ^ c between <5iV c and SN, we find that the term in SN C is 
suppressed compared to the corresponding term in SN by 

-^(Vxx ~ Vvv) x (Vvv or 77 xx )x(sin 2 7 or cos 2 7). (43) 

Here — ip 2 + x 2 . Since we are working in the vacuum dominated regime and slow-roll, Eq. ([3]), p-?7 2 v <C 1 and 

2 9 1 1 
t-J/vv ^ 1- We therefore see that SN C is greatly suppressed compared to SN. 

9i xx 

D. Further evolution after inflation 

So far in this section we have assumed a quick transition to the radiation epoch at the end of inflation, thereby 
neglecting the dynamics of the waterfall field. However if we consider the role of the waterfall field, then after the 
waterfall field is destabilised there may be a further evolution of the primordial curvature perturbation, which will 
lead to a change of the observable parameters. This applies to any model with an inhomogencous end of inflation 
since there are isocurvature perturbations still present after the waterfall field is destabilised and inflation has ended. 
Further evolution will depend on the details of reheating in a model dependent way. To the best of our knowledge 
this issue has not been considered in depth in any paper on an inhomogeneous end of inflation. If we assume an 
instantaneous transition to radiation domination (so a completely efficient and immediate decay of the waterfall and 
inflaton fields) then there will be no further change to the observables as we have argued in the previous section. 
However this is clearly an idealised case. For a review of reheating after inflation see for example [49j . 

In the special case where the waterfall field is also light during inflation Barnaby and Cline [3| have shown there is 
the possibility of generating a large non-Gaussianity during preheating for certain parameter values. This is possible 
even if there is only one inflaton field and the waterfall field present. However in this case inflation does not end 
abruptly when the waterfall field is destabilised so this is not the scenario we have considered in this paper. 

In practice the efficiency of reheating will depend on the couplings between the waterfall and inflaton fields to 
any preheat fields, as well as on the ratio between 51,(72 and A [50j. In one regime where A 3> gi and g2 preheating 
depends mainly on the inflaton fields. In the case where the two inflaton fields couple identically to all further particles 
preheating depends on the coupling constants g\ and (72 and only much more weakly on their bare masses proportional 
to r\ vv and ?7 XX j50l |. So in the case where g\ = f/2 which we have considered earlier it may be reasonable to expect 
little or no further evolution of the curvature perturbation, because the isocurvature perturbations should be irrelevant 
and decay during preheating. On the other hand if g\jgi <C 1 which is the case considered in some previous works 
on an inhomogeneous end to inflation the effect from preheating is more likely to be important. In another regime 
where A -c g\ and g 2 preheating depends mainly on the waterfall field, so it may be that in this case the isocurvature 
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perturbations are again unimportant. However we should also consider the period from when the waterfall field is 
destabilised and until preheating begins. It could be that the amount of expansion during this time also leads to an 
extra correction to the curvature perturbation. This is a complicated issue which deserves further investigation but 
is beyond the scope of this paper. To the best of our knowledge there has been no study of preheating or reheating 
in a model of hybrid inflation with more than one inflaton field. 



IV. SCALE DEPENDENCE 



In general the bispectrum parameterised by /nl can have both a scale and a shape dependence. We are considering 
the local form of /nl which means that /nl is shape independent. However it can still have a (slow- roll suppressed) 
scale dependence [U HH EH . 

In our examples /nl has a scale dependence both because of the exponential term in /nl, (|17p . and because f will 
vary through the change of the initial value of sin 2 0* . We find 



dink dink 
Using this we find from p7|) that 



-2(r] vv - i lxx ). (44) 



n /NL - i = = _ 4 2ap3a. (45) 

/NL d log k 1 + r y ' 

In the case that we include the effect from the surface where the waterfall field is destabilised and g\ = g\ we find 
from ([M]) that 

»*l - 1 = -4 ^Zpm . (46 ) 

For both cases the spectral index of /nl satisfies 

- tijbw ~Vxx)< "/nl - K 0, (47) 

for any value of r and hence /nl will be smaller on small scales as we can see in the figures [T] and [5J 

Because we require a relatively large value of r\ w — r\ xx > 1/N for our model to generate a large non-Gaussianity 

it is quite possible for our model to generate a relatively significant scale dependence of /nl- However the amount 

also depends on f and when this is large then n / NL — 1 is suppressed. 

We note that this is in contrast to the large non-Gaussianity from an inhomogeneous end of inflation found in [2l| . 

In the specific cases they considered to generate a large non-Gaussianity the non-Gaussianity was generated purely 

at the end of inflation and /nl is scale independent. In detail we see from Eqs. (4.4) and (4.24) in [2l[ that their 

formulae for /nl does not depend on N or on any quantities evaluated at Hubble exit. 



V. CONCLUSION 



We have made an in depth study of a model of hybrid inflation with two inflaton fields (two-brid inflation) . We 
have studied the parameter space where non-Gaussianity may be large during slow-roll inflation. In particular we 
have calculated the observable parameters Vq, tlq — 1, r, /nl, 9nl and tnl- The spectral index depends on the 
weighted sum of the slow-roll parameters r) vv and r\ xx while the bispectrum depends exponentially on the difference 
of the same slow-roll parameters, so it is possible to have a very large bispectrum and a scale invariant spectrum. We 
have shown that the trispectrum is generally large through tjvl whenever the bispectrum is large, but that the other 
parameter which parameterises the trispectrum, g^h is always smaller than /nl and strongly suppressed compared 
to tnl- We have also shown that for certain initial conditions it is possible that tjvl is the only large non-linearity 
parameter, so the first observational signature of non-Gaussianity from this model could be the trispectrum through 
tnl- 

Furthermore we have shown that during slow-roll inflation the loop corrections to the power spectrum, /nl and 
are always suppressed compared to the tree level terms. The suppression follows from the constraint on the 
background trajectory that the classical background value of the field should dominate over its quantum fluctuations. 
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We then investigated how the large non-Gaussianity which was generated during slow-roll may be changed by the 
effects from the end of inflation. We included the waterfall field which ends inflation when its effective mass becomes 
negative and it is destabilised. The effect of the waterfall field depends on the values of the coupling constants between 
the waterfall field and the two inflaton fields. If we choose the coupling constants so that surface where the waterfall 
field is destabilised corresponds to a hypersurface of uniform energy density then there is no change to the observables 
we calculated during slow roll. However if we choose the two coupling constants to be equal then observables at both 
linear and higher order arc changed by an amount that depends on the ratio r) vv /i] xx . 
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VI. APPENDIX 



A. Loop Suppression 

We require that the quantum fluctuation does not overwhelm the classical evolution of the fields. In order to do 
that we require \x\tt/H 2 > y/3/2 for both fields on the background trajectory With this constraint we show 

that the tree level dominates the loop correction. In the slow-roll condition, this constraint leads to 



/ / 



Tjipip 



/ / 



Vxx 



(48) 



or equivclantly 



< 



\/6 



$X* 



< 



Vxx 



V6 



(49) 



Using the 6N formalism, 



(50) 



in hybrid inflation, we find the first derivatives of the number of e-foldings 



-2Nri vv 



2r 



Y v 



-2JVr, xx 



2e 



(51) 



and the second derivatives 



n =*k 



N. 



-4?7 



xx 



2 l ] 



(52) 



vx 



2N v N x [ p - (Wipv+Vxx) 



where 



7 _ rftyiple-™"™ +y xx xle- 2Nr >™ 



(53) 
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Wc note that 

p <Vmax = max{r) lpv ,r) xx }. (54) 
From the definition of the power spectrum and observations 
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(55) 



we know that 

\N v S(p\ ^ \N V H*\ < 1(T 4 , \N x 5 X \^\N x H*\<10- i . (56) 

Using Eq. (|4"9"|) and (|56"|) we can easily check that the second order term in the expansion of SN are always smaller 
than the first order term. For example 



l -N w {5^) 2 = \^{5^) 2 + f-2 % + J] (N V S^) 2 < I (^) (N V S^) + r, max lO- 5 {N v S^) < \r, vv N v \S^, 

Z 2 If* \ p J A \ Z7T / 

(57) 



which is similar for the N xx term and 

N vx (SipJ X *) = 2{N v 5<p*)(N x 5 X *) Q " (v vv + Vxx)j < 10- 4 |?/ ma:c |max(|iV v |^, \N X \6 X *). (58) 

The possible one loop domination can happen when \ 1S much smaller than S\*, and (r? V¥ , — r] xx )N > 1 which is 
the case of Yeinzon ct al. [22[ , but this can be possible only when we break the classical constraint we supposed in 
Eq. jp). 

In \22\, [23[ the dominant loop correction to the power spectrum, bispectrum and trispectrum was found to come 
from a single loop diagram in each case which is the loop expansion one finds if truncating the SN expansion at second 
order. Specifically the loop corrections they found which may be dominant are (see figure H]) 



Pl l °° P N AB N AB 

ptree (N C N C ) 2 



V (59) 



ri loop 5 N AB N BC N AC 

f ™ = 6 (N D N D f P <> (60) 

Hoop N ab N bc N C dN ad 

TNL = (N E N E f (61) 

Similar formula for /nl and tnl at tree level are given in (|19| (figure [3|). By comparing the terms at tree and 1-loop 
level (which are suppressed by a factor of ~ 10 -10 ) it is clear that the loop terms can only be large in the case 
that one of the second derivative terms in the SN expansion is extremely large. The integral over the loop momentum 
gives rise to a logarythmic infrared divergence ln(fcL) where L is the large scale cut off. We have followed [12, H3, 5(| 
in assuming that we can take \n(kL) ~ 1. In the terminology of 15411 this corresponds to working in a minimal box. 
For more discussion of this divergence see for example |45l l55l . l56l. l57j . 

Without loss of generality we assume \N vip \ < max(|N vx |, |N XX |), so we need at least one of {N^l/N 2 > 1 or 
\N XX \/N 2 ^> 1 in order to have a significant loop correction. First we see that the cross derivative term is never large, 
using Eq. (J3IJ) and (|52|) . 

I AT I ,/f 

1 vxl <i-^-_\r, max \<2\r, max \. (62) 



N$ + Nl 1 + r 

The ratio which can be large is N xx /N 2 , in the case that it is large we have the simplified formulas 



P l l ° OP „ (N xx \ 2 1 



ptree ' \ AT2 
C \ 1 V 



tree ' \ AJ2 I /-i , ~\2 



(1 + r) 



z^Pc, (63) 
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j-tree ^ 5 N xx f j loop ^ 5 / -/V^A 1 ^ /^..n 

/m - 6 JV| (l+jf' , » L -6\N}) (TTff 

We recall the definition r = (N x /N ip ) 2 . So it appears that by making N xx /N 2 large enough we can make the loop 
corrections as large as we like. 

However using Eq. (|56p and (|57p . we find the inequality 



NxA (66) 



2 

n 2 y f(i + f) " '' xx 

and hence we see from (|63[) f|64[) and (|65|) that the loop correction to the power spectrum, /nl and tjvl is always 
suppressed. 

B. Other loop corrections 

So far we have proved that the loop corrections which [22], HH found to be largest are still suppressed compared to 
the tree level terms. Here we show that these particular loop corrections remain the dominant loop corrections even 
for a more general trajectory in field space. To do this we need to go beyond second order in the SN expansion. In 
fact to prove that all of the 1-loop corrections to the power spectrum, /nl and t^l are suppressed we need to go up 
to fourth order in the SN expansion. 

At third order, the only derivatives which can be very large are 

Ntpxx - 2n ^* Nxxx - 6s/p(n - 2n )^^- (67) 

J V ip lyl if iV lp if 

while at fourth order the largest derivative is given by 

^ = 6(, w -2, xx )(^) 2 . (68) 

We use the third order derivatives to calculate <7tvl, see (TTpl) and (|2"Tj) . 

For details of the loop corrections we refer the reader to [45| , in which a diaggrammatic approach to calculating the 
primordial n-point functions of £ at any required loop level was developed. The formula for the n-point functions 
depend on the derivatives in the SN formalism, for every vertex with r internal legs there is a corresponding term 
with r derivatives of N. We note from Fig. [3] that the tree level diagrams for the power spectrum, /nl and tml 
all depend on diagrams with vertices connected to either one or two internal lines. The dominant loop correction 
found by (22I |23| in each case follows by adding an internal line between the two vertices with a single internal line 
attached in the tree level diagrams. These diagrams are shown in Fig. 0] This converts two first order derivatives of 
N into second order derivatives, which adds a large multiplicative factor of (N xx /N 2 ) 2 to the 1-loop term. When 
constructing any other 1-loop diagram we can either dress a vertex, which means adding two internal lines to a single 
vertex or connect an internal line between two external vertices where at least one of them does not have a single 
internal line atteched at tree level. In each of these cases it follows from (f6"T)) and (f6"8"]) that these other diagrams are 
not boosted by a factor larger than N xx /N 2 and hence are subdominant to the loops found in [22l.[23|. 

We conclude that all one-loop terms to the power spectrum, /nl and tjvl are suppressed compared to the tree 
level terms, even for a general trajectory in field space. In fact this argument applies to any loop term, provided that 
we can truncate the expansion of 5N at fourth order. When drawing an I + 1-loop level diagram we add an extra 
internal line to a diagram at Z-loop level. This means attaching in total two extra derivatives to the derivatives of 
N corresponding to the vertices where the extra internal line is added. At most this can add a numerical factor of 
(N xx /N 2 ) 2 . However there is also a suppression factor of order Vq which comes from going to a higher order in loops, 
which as we have seen explicitly in the case of the one loops diagrams always leads to a suppression. 
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FIG. 3: Tree level diagrams for the power spectrum (left hand side), the bispectrum (centre) and the relevant tree diagram for 
the trispectrum (right) which corresponds to tnl- The diagrams were drawn using JAXODRAW [k| 



'I ll 

I I I 

I l I 

-J I J_ 



FIG. 4: The dominant one-loop level diagrams for the power spectrum (left hand side), the bispectrum (centre) and the 
trispectrum (right). 
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